Search for: All records

Abstract Collapsars—rapidly rotating stellar cores that form black holes—can power gamma-ray bursts and are proposed to be key contributors to the production of heavy elements in the Universe via the rapid neutron capture process (r-process). Previous neutrino-transport collapsar simulations have been unable to unbind neutron-rich material from the disk. However, these simulations have not included sufficiently strong magnetic fields and the black hole (BH), both of which are essential for launching mass outflows. We presentνh-amr, a novel neutrino-transport general relativistic magnetohydrodynamic (νGRMHD) code, which we use to perform the first 3D globalνGRMHD collapsar simulations. We find a self-consistent formation of a weakly magnetized dense accretion disk, which has sufficient time to neutronize. Eventually, substantial magnetic flux accumulates near the BH, becomes dynamically important, leads to a magnetically arrested disk (MAD), and unbinds some of the neutron-rich material. However, the strong flux also hinders accretion, lowers density, and increases neutrino-cooling timescale, which prevents further disk neutronization. Typical collapsar progenitors with mass accretion rates, $\dot{M}\sim 0.1-1M_{\odot }{\mathrm{s}}^{-1}$, do not produce significant neutron-rich (Y_e < 0.25) ejecta. However, we find that MADs at higher mass accretion rates, $\dot{M}\gtrsim \text{few}{M}_{\odot }{\mathrm{s}}^{-1}$(e.g., for more centrally concentrated progenitors), can unbindM_ej ≲ M_⊙of neutron-rich ejecta. The outflows inflate a shocked cocoon that mixes with the infalling neutron-poor stellar gas and raises the final outflowY_e; however, the finalr-process yield may be determined earlier at the point of neutron capture freeze-out. Future work will explore under what conditions more typical collapsar engines becomer-process factories. 

Abstract We introduceMahakala, aPython-based, modular, radiative ray-tracing code for curved spacetimes. We employ Google’sJAXframework for accelerated automatic differentiation, which can efficiently compute Christoffel symbols directly from the metric, allowing the user to easily and quickly simulate photon trajectories through non-Kerr spacetimes.JAXalso enablesMahakalato run in parallel on both CPUs and GPUs.Mahakalanatively uses the Cartesian Kerr–Schild coordinate system, which avoids numerical issues caused by the pole in spherical coordinate systems. We demonstrateMahakala’s capabilities by simulating 1.3 mm wavelength images (the wavelength of Event Horizon Telescope observations) of general relativistic magnetohydrodynamic simulations of low-accretion rate supermassive black holes. The modular nature ofMahakalaallows us to quantitatively explore how different regions of the flow influence different image features. We show that most of the emission seen in 1.3 mm images originates close to the black hole and peaks near the photon orbit. We also quantify the relative contribution of the disk, forward jet, and counterjet to 1.3 mm images. 

Abstract Long-duration gamma-ray bursts (lGRBs) originate in relativistic collimated outflows—jets—that drill their way out of collapsing massive stars. Accurately modeling this process requires realistic stellar profiles for the jets to propagate through and break out of. Most previous studies have used simple power laws or pre-collapse models for massive stars. However, the relevant stellar profile for lGRB models is in fact that of a star after its core has collapsed to form a compact object. To self-consistently compute such a stellar profile, we use the open-source code GR1D to simulate the core-collapse process for a suite of low-metallicity rotating massive stellar progenitors that have undergone chemically homogeneous evolution. Our models span a range of zero-age main-sequence (ZAMS) masses:M_ZAMS= 13, 18, 21, 25, 35, 40, and 45M_☉. All of these models, at the onset of core-collapse, feature steep density profiles,ρ∝r^−α, withα≈ 2.5, which would result in jets that are inconsistent with lGRB observables. We follow the collapses of four of the seven models until they form black holes (BHs) and the other three models until they form proto-neutron stars (PNSs). We find, across all models, that the density profile outside the newly formed BH or PNS is well represented by a flatter power law withα≈ 1.35–1.55. Such flat density profiles are conducive to the successful formation and breakout of BH-powered jets and are, in fact, required to reproduce observable properties of lGRBs. Future models of lGRBs should be initialized with shallower post-collapse stellar profiles, like those presented here, instead of the much steeper pre-collapse profiles that are typically used. 

Abstract We combine our dynamical modeling black-hole mass measurements from the Lick AGN Monitoring Project 2016 sample with measured cross-correlation time lags and line widths to recover individual scale factors,f, used in traditional reverberation-mapping analyses. We extend our sample by including prior results from Code for AGN Reverberation and Modeling of Emission Lines (caramel) studies that have utilized our methods. Aiming to improve the precision of black-hole mass estimates, as well as uncover any regularities in the behavior of the broad-line region (BLR), we search for correlations betweenfand other AGN/BLR parameters. We find (i) evidence for a correlation between the virial coefficient ${\log }_{10}\left(f_{\mathrm{mean},\sigma }\right)$and black-hole mass, (ii) marginal evidence for a similar correlation between ${\log }_{10}\left(f_{\mathrm{rms},\sigma }\right)$and black-hole mass, (iii) marginal evidence for an anticorrelation of BLR disk thickness with ${\log }_{10}\left(f_{\mathrm{mean},\mathrm{FWHM}}\right)$and ${\log }_{10}\left(f_{\mathrm{rms},\mathrm{FWHM}}\right)$, and (iv) marginal evidence for an anticorrelation of inclination angle with ${\log }_{10}\left(f_{\mathrm{mean},\mathrm{FWHM}}\right)$, ${\log }_{10}\left(f_{\mathrm{rms},\sigma }\right)$, and ${\log }_{10}\left(f_{\mathrm{mean},\sigma }\right)$. Last, we find marginal evidence for a correlation between line-profile shape, when using the root-mean-square spectrum, ${\log }_{10}{\left(\mathrm{FWHM}/\sigma \right)}_{\mathrm{rms}}$, and the virial coefficient, ${\log }_{10}\left(f_{\mathrm{rms},\sigma }\right)$, and investigate how BLR properties might be related to line-profile shape usingcaramelmodels. 

Abstract We have modeled the velocity-resolved reverberation response of the Hβbroad emission line in nine Seyfert 1 galaxies from the Lick Active Galactic Nucleus (AGN) Monitoring Project 2016 sample, drawing inferences on the geometry and structure of the low-ionization broad-line region (BLR) and the mass of the central supermassive black hole. Overall, we find that the HβBLR is generally a thick disk viewed at low to moderate inclination angles. We combine our sample with prior studies and investigate line-profile shape dependence, such as ${\log }_{10}\left(\mathrm{FWHM}/\sigma \right)$, on BLR structure and kinematics and search for any BLR luminosity-dependent trends. We find marginal evidence for an anticorrelation between the profile shape of the broad Hβemission line and the Eddington ratio, when using the rms spectrum. However, we do not find any luminosity-dependent trends, and conclude that AGNs have diverse BLR structure and kinematics, consistent with the hypothesis of transient AGN/BLR conditions rather than systematic trends. 

Abstract We carried out spectroscopic monitoring of 21 low-redshift Seyfert 1 galaxies using the Kast double spectrograph on the 3 m Shane telescope at Lick Observatory from 2016 April to 2017 May. Targeting active galactic nuclei (AGNs) with luminosities ofλL_λ(5100 Å) ≈ 10⁴⁴erg s⁻¹and predicted Hβlags of ∼20–30 days or black hole masses of 10⁷–10^8.5M_⊙, our campaign probes luminosity-dependent trends in broad-line region (BLR) structure and dynamics as well as to improve calibrations for single-epoch estimates of quasar black hole masses. Here we present the first results from the campaign, including Hβemission-line light curves, integrated Hβlag times (8–30 days) measured againstV-band continuum light curves, velocity-resolved reverberation lags, line widths of the broad Hβcomponents, and virial black hole mass estimates (10^7.1–10^8.1M_⊙). Our results add significantly to the number of existing velocity-resolved lag measurements and reveal a diversity of BLR gas kinematics at moderately high AGN luminosities. AGN continuum luminosity appears not to be correlated with the type of kinematics that its BLR gas may exhibit. Follow-up direct modeling of this data set will elucidate the detailed kinematics and provide robust dynamical black hole masses for several objects in this sample.